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^ Abstract. It is well established that phenomenological two-zone diffusion models 

of the galactic halo can very well reproduce cosmic-ray nuclear data and the observed 
qq antiproton flux. Here, we consider lepton propagation in such models and compute the 

expected galactic population of electrons, as well as the diffuse synchrotron emission 
\q that results from their interaction with galactic magnetic fields. We find models in 

<^> agreement not only with cosmic ray data but also with radio surveys at essentially all 

frequencies. Requiring such a globally consistent description strongly disfavors very 
. . large (L > 15 kpc) and, even stronger, small (L < 1 kpc) effective diffusive halo sizes. 

^ This has profound implications for, e.g., indirect dark matter searches. 

'X 



1. Introduction 



A wealth of observational data strongly suggests that diffusion governs the propagation 
of galactic cosmic rays (CRs) pQ. Any model for the underlying physical processes 
requires a basic assumption about the geometry of the region the CRs sweep. In diffusive 
models, the Galaxy is represented by a thin disk sandwiched by a thick magnetic diffusive 
halo with cylindrical symmetry. Given the intrinsic limitations to any such model, as 
well as present-day CR data, a considerable uncertainty in the propagation parameters 
is generally unavoidable. 

A realistic 3D modeling turns to extensive computer codes, such as Galprop [21 [3] or 
Dragon [I], aiming at a detailed description of sources, interstellar matter distribution, 
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magnetic field structure and diffusion phenomena. While reaching a high degree of 
accuracy - particularly needed for the gamma-ray component - such an approach is 
computationally expensive and does not always make it straight-forward to extract 
physical answers for the underlying processes and input parameters. An effective 2D 
approach, on the other hand, benefits from analytical solutions to the spatial diffusion 
equation [3 El [7j, which allows fast computations and thus efficient scans of the 
parameter space [8]. Despite a small number of free parameters, it can consistently 
describe both nuclear CRs [5], CR antiproton [9] as well as electron and positron data 



Here, we investigate whether even radio data can be understood in this 
2D framework. We start by inferring the galactic electron distribution from an 
interpretation of large-scale radio survey data as synchrotron radiation and then 
compare this to the expectation in our propagation model. We find that radio 
observations are indeed consistent with models that correctly describe CR data, 
providing a remarkable hint that our effective propagation model is not too far from a 
real picture of galactic phenomena, at least on kpc scales. Furthermore, we show that our 
procedure can be used , in principle, as a new method to constrain propagation models 
which is complementary to using the boron over carbon (B/C) ratio [5] or radioactive 
isotopes ratios (i.e. 10 Be/ 9 Be) 00 E]. 

The structure of this article is as follows. We start by reviewing, in Section [2j 
how relativistic electrons produce synchrotron radiation when propagating through the 
galaxy. In Section [3j we compute the galactic electron distribution in our diffusion 
model and compare the expected synchrotron radiation in Section [4] to radio surveys 
at various frequencies, demonstrating that synchrotron radiation indeed provides a very 
promising means of both inferring properties of the interstellar electron distribution and 
to provide constraints on the adopted diffusion model. After a discussion of possible 
biases in our analysis in Section [5j we present our conclusions in Section [6j 

2. Synchrotron radiation 

Relativistic electron^] emit synchrotron radiation while propagating through the galactic 
magnetic field [T2] . For electrons with energy E and a magnetic field of strength B, the 
emission power per unit frequency is given by 



where v c = 3eBE 2 / (inmlc 5 ) , F(x) = x f^°d( K 5 / 3 (() and K 5 / 3 is a modified Bessel 

— * 

function; we take the average over the angle 9 between the electron momentum and B 
because we assume an isotropic electron distribution (note also that, in general, not only 
the regular but also the turbulent component of the galactic magnetic field contributes 
to the total signal - which further motivates this average). 

| Unless explicitly stated otherwise, we will in the following use the term electron to denote both 
electrons and positrons. 



noun]. 




(i) 
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For a given electron density n e (E,x), the expected intensity in synchrotron 
radiation is thus given by 



where J u = J ^ ^ dE is the emissivity, a u the inverse of the absorption length and 
the integrations are taken along the line of sight towards the observed direction in 
the sky, averaged over an angular region AQ. Since radio emission is often associated 
with thermal phenomena, the intensity is traditionally stated in terms of temperature: 
with the Rayleigh- Jeans law in mind, I v = 2u 2 k B T/c 2 , one can define the brightness 
temperature as 



The galactic magnetic field has an average strength of G(/j,G) [13]. While its 
detailed structure in reality can be expected to be rather sophisticated [H] , we will here 
adopt an effective approach and treat it to be spatially constant within the diffusion zone 
as we will only be interested in an average, large-scale description of radio data above 
the galactic plane. In fact, this simplifying assumption seems necessary in order to be 
consistent with the homogeneous diffusion coefficient that enters as a basic ingredient 
to our propagation model. We verified that adopting instead a magnetic field falling off 
exponentially away from the disk - in principle also consistent with the geometry of our 
diffusion model - does not change significantly our predictions for the synchrotron flux 
(integrated along the line of sight) with respect to our choice of a constant magnetic 
field that abruptly vanishes at the border of the diffusion zone. This small difference 
can readily be understood in terms of the propagation scale length, which is ~ 300 pc 
for electrons with energies ~ 3 GeV: most of the electrons will simply remain relatively 
close to the galactic disk, where the magnetic field distributions are similar, thereby 
reducing the differences in the integrated synchrotron emission. 

For our analysis, we take into account absorption of synchrotron photons by 
both thermal electrons and ions from a possible hot gas component in the disk 
[15J. Synchrotron self-absorption is another potential effect, but for B ~ fiG it 
becomes numerically important only for frequencies below 1 MHz [TU] where it is 
however dominated by free-free absorption. Finally, let us mention that for an electron 
distribution following a power law, dn e /dE oc E~ J , Eqs. ([l](3]) tell us that also the 
synchrotron intensity, in the case of negligible absorption, follows a power law, T b oc u~ a , 
with a spectral index 




(2) 




(3) 



a = 



(7 + 3)/2. 



(4) 



3. Galactic electron population 



The transport parameters of the two-zone diffusion model introduced in Ref. [5] are 
determined from the B/C analysis and correspond to the size of the diffusive halo of the 
Galaxy L, the normalization of the diffusion coefficient K and its slope 5 (defined by 



Radio data and synchrotron emission in consistent cosmic ray models 



4 



Model 


prop 

L [kpc] 


. paramet( 

A o Myr 


3rs 

5 


radio dat 

408 MHz 


a (x 2 /cl.o.f.) 
1.42 GHz 


min 


1 


0.0016 


0.85 


11.6 (6.8) 


11.9 (6.3) 


med 


4 


0.0112 


0.70 


4.9 (2.0) 


4.9 (2.0) 


max 


15 


0.0765 


0.46 


10.8 (4.8) 


8.9 (3.9) 



Table 1. Benchmark models compatible with B/C data [17]. Both 'min' and 'max' 
are clearly disfavored by radio data towards the galactic anti-center (averaged over 10° 
(15°) and excluding the disk). 

K = K (3R S , where R = p/q is the rigidity), with a rather strong degeneracy between 
the allowed parameters - in particular between Kq and L. Diffusive reacceleration and 
convection, while crucial for the nuclei analysis, have been shown to be only mildly 
relevant for lepton fluxes [in]. Indeed, the two processes shape GeV electron fluxes in 
the opposite direction, so that their combined effect is rather small; for that reason, 
and for the sake of simplicity, we do not take them into account here. In this seminal 
analysis, we will for simplicity mainly refer to the three benchmark propagation models 
shown in Table [TJ among all models compatible with B/C data, these were shown to give 
the minimal, medium and maximal flux in antiprotons, respectively, that is expected 
from dark matter annihilations in the galactic halo [T7] . 

High energy electrons are produced in galactic accelerators such as supernova 
remnants (SNRs) or pulsars (primary electrons), as well as in hadronic interactions of 
galactic protons and helium nuclei with the interstellar medium (secondary electrons). 
We calculate the primary e~ flux from SNRs following Ref. [TT] and the subdominant 
secondary e ± component as described in Ref. [TO], using a cylindrical gas density 
distribution in the disk as given in Ref. [H] . For the energy losses, we take into account 
inverse Compton scattering off the interstellar radiation field, as well as synchrotron, 
bremsstrahlung and ionization losses in the interstellar medium [8]. 

4. Comparison to radio data 

Starting from the 1960s, several groups have performed large-area radio surveys in the 
frequency range from about 1 MHz to 100 GHz; for an extensive list, we refer the 
reader to Ref. [Ej. In this reference, the surveys at 0.010 [20J, 0.022 [21], 0.045 [22], 
0.408 [23], 1.42 [21], 2.326 GHZ [25], as well as WMAP foregrounds at 23, 33, 41, 61 
and 94 GHz [20], were transformed to galactic coordinates and pixelized. In Fig. [TJ we 
show the north and south pole at these frequencies, with the isotropic CMB component 
removed and averaged over a circular region with 10° in diameter. 

At these frequencies and latitudes, the dominant source of the radio signals should 
be synchrotron radiation (except for the excess seen in the 94 GHz band which is 
probably due to spinning dust [27] and which we will not consider in the following). Free- 
free emission in general becomes important at frequencies v > 1 GHz, but is physically 
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Figure 1. Observed radio fluxes vs. frequency. Data points are generated using the 
GSM software [T3] which exactly reproduces the observational data at 10 MHz [20] . 
22 MHz [2T], 45 MHz [22], 408 MHz [2J, 1.42 GHz [24 , 2.326 GHz [25], as well as 
WMAP foregrounds at 23, 33, 41, 61 and 94 GHz [25]. The thick solid line represents 
the best fit for synchrotron radiation from an electron population (after propagation) 
with dn e /dE oc E~ lx (dn e /dE oc £T 72 ) for E < E hl (E > E hT ); the dashed lines 
indicate the effect of changing the best-fit value E^ T = 3.9 GeV. The thin solid line 
shows a fit with a break in the energy losses for the 'med' model of Tab. [I] see text 
for further details. 

subdominant when looking away from the galactic disk: adopting values for the thermal 
electron temperature and distribution as given in Ref . [28] , we estimate the highest free- 
free contribution (at 61 GHz) to be less than 20% in a cone towards the poles; given the 
error bars shown in Fig. [TJ this does not have a significant effect on our analysis. From 
the data, one can clearly distinguish two regimes with a different power-law behavior 
in frequency - which directly translates to the necessity of a break in the spectral 
index of the propagated galactic electron population. If the propagated electrons are 
simply modeled with a broken power law, the best-fit values for dn e /dE are a break 
in the spectral index A7 = 1.57^25 at ^br = 3.87±^ GeV, with 72 = 3.45j£aj above 
E bT (note that the sharp break in the electrons gets smoothed because we use the full 
expression ([!]) for the synchrotron power rather than the often adopted monochromatic 
approximation). This result is fully consistent with, e.g., that of Strong et al. [36J, see 
their Fig. 5, and displayed as a thick, black solid line. Here, we chose a fiducial value of 
B = 6.5 fiG; a different value would simply change the location of E br and the (arbitrary) 
overall normalization (by a factor oc B 2 ), but not the functional dependence on v (we 
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Figure 2. Measured spectrum in CR electrons (e~ for the case of PAMELA) below 
50 GeV [521 [301 131] and predictions according to Refs. (TUJ [TT], assuming a spectral 
break in the energy losses. This spectrum is consistent with the radio data shown in 
Fig. [I] (vertical bands indicate the corresponding frequency of synchrotron radiation) . 

take the opportunity to remind the reader that B = 6.5 fiG is an effective value which 
takes into account both the average regular and turbulent component of the galactic 
magnetic field). 

Such an electron distribution is actually also in rather good agreement with the 
observed CR electron fluxes (note that only electrons well below 100 GeV are relevant 
to our discussion because the spectrum at higher energies is likely to be dominated 
by local sources). In Fig. [2j we show the Pamela e~ data [29J together with the 
e = e + + e~ data taken by Fermi-LAT |30j . Very recently, the Fermi-LAT Collaboration 
has performed a separation between observed positrons and electrons above 20 GeV 
using the Earth magnetic field as a spectrometer [31 J, confirming their previous results 
on the total lepton fluxes. 

Along with the data, we plot in Fig. [2] theoretical predictions for the interstellar 
and solar modulated fluxes. Black lines correspond to the radio best fit (see Fig. [TJ, 
where electrons are naively shaped by two power laws with break of A7 = 1.57. We 
also show the resulting flux at Earth after propagation according to Refs. pm EH] and 
for the propagation models listed in Table [T] here, we included for comparison a break 
of 1.1 in the energy loss term, making the hypothesis that losses effectively follow E 0,9 
below Ebreak (for the 'med' model, this is in addition shown as a red solid line in Fig. [TJ. 
In order to obtain the curves shown in Fig. [2j we choose the normalization and spectral 
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index 7i n j of the primary e~ spectrum in such a way as to fit the total (primary and 
secondary) e~ spectrum to the low-energy PAMELA data (we checked that adding the 
Fermi electron data, given their energy range and error bars, would not significantly 
affect this normalization). We stress that our normalization is consistent with Fig. [TJ 
i.e. it reproduces both the the radio data (with B = 6.5 /iG as stated above) and, 
roughly, the low-energy lepton data [29] . 

For energy losses that do not significantly change during the typical path an electron 
propagates, the spectral index of the expected electron flux is approimately given by 

m 

7«7inj + ^ + 5-l), (5) 

where 7i n j is the spectral index of the injected electrons (before propagation) and the 
energy loss term is assumed to scale like b = —d t E oc E 13 . This relation is not more 
valid, however, if electrons experience different energy loss regimes while wandering 
through the Galaxy, depending for example on the time spent in the halo (which is 
gas free) or in the disk. If energy losses are dominated by a term specific to the local 
environment where the electron spends most of its time, b oc E^ oc , one can derive from 
the general expressions given in Ref. [UJ that the expected relation for the observed 
spectrum instead becomes 

7~7inj + ^(-/3 + 5-l) + Aoc, (6) 

which reduces to Eq. ^ for (3\ oc = (3. We note therefore that the break illustrated in 
Figs. [TJ and [2j taking into account the viable scatter in A7, might simply reflect the fact 
that at lower energies the scattering on thermal ions and electrons (6] oss oc E) becomes 
more important than inverse Compton losses (&i OS s oc E 2 ). This interpretation is also 
consistent with the location of the break that is determined by the relative strength of 
these processes: taking into account uncertainties in the radiation, magnetic field and 
gas densities, we expect the transition to occur in the range 1 GeV< inbreak < 10 GeV. 
In fact, as also shown in Figs. [TJand[2j a corresponding value of A(3\ oc = 1.1 (which we 
choose slightly larger than 1 in order to capture possible effects of reacceleration) does 
show a very reasonable agreement with both radio and lepton data. 

Let us now turn to the angular shape of the radio signals, which we show in Fig. [3] 
for the case of the Haslam [23] map at 408 MHz. In the same figure, we indicate the 
expectation for the models of Tab. [TJ including the full variation of the parameters K 
and 5 within the range compatible with the B/C analysis performed in Ref. [5]. One can 
clearly see that the size of the diffusive halo has a rather strong impact on the angular 
shape of the resulting synchrotron signal. Demanding consistency with B/C data, small 
halo sizes L ~ 1 kpc are essentially excluded, but also large values L > 15 kpc show some 
tension with radio data. Note that the magnetic field normalization does not affect the 
angular shape of the synchrotron emission; here, it was chosen such as to be consistent 
with Figs. [TJ and [2j 
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Figure 3. Brightness temperature versus latitude at 408 MHz. Solid lines correspond 
to the propagation models given in Table [T] The shaded areas are obtained by varying 
Kq and 5, at fixed L, over the full range compatible with B/C data. 



In order to demonstrate the constraining power of the angular shape of radio data 
for diffusion models, we next treated the overall normalization of the synchrotron signal 
(aka the magnetic field) as a free parameter to minimize \ 2 f° r the min/med/max 
models. The result is shown in Tab. [T] for the 408 MHz and 1.42 GHz maps; as 
anticipated, it is very difficult to reconcile radio data with the max and, even more 
so, the min model. Similar conclusions apply to all surveys with v < 408 MHz, while for 
v > 1.42 GHz the data generally start to reproduce the expected synchrotron pattern 
worse, even for the med model; this is likely due to molecular clouds and pulsars (and/or 
SNRs), mostly in the southern hemisphere, that produce bremsstrahlung emission. We 
note that even the overall fit quality, for the min model, is quite good for angular averages 
of at least 15°; on smaller scales, on the other hand, we can obviously not expect our 
effective model to reproduce the detailed features that are visible in the radio sky. While 
these results already indicate the power of the method and thus warrant a more detailed 
analysis of the propagation parameter space and the magnetic field structure [32], this 
is beyond the scope of the present work. 

5. Discussion 

Let us now discuss which assumptions in the analysis presented above are affected 
by uncertainties that might bias our conclusions. We have explicitly verified that the 
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absorption of synchrotron radiation by thermal electrons is negligible for v > 10 MHz 
and insensitive to the thermal electron temperature (for 4000 K< T e <8000 K) below this 
frequency; for very low values (T e <3000 K) absorption modifies the spectral shape in a 
way strongly disfavored by the data. Weak uncertainties are derived from the ensuing 
bremsstrahlung radiation at high frequencies, whose contribution emerges around the 
GHz and gets stronger with lower T e . 

We checked that a break in the diffusion coefficient does not induce relevant changes 
in the frequency spectrum, in contrast to the break in the (total) propagated spectrum 
discussed above; we found that a break in the source spectrum of primary electrons 
alone, on the other hand, would result in an overproduction of synchrotron radiation 
from secondary electrons at frequencies v < 100 MHz. We also verified that modifying 
K or L merely shifts the spectrum w.r.t. v by changing its normalization; as a function 
of u, the three benchmark models of Tab. [TJ e.g., simply scale up the spectrum by a 
factor of 5 when going from min to max. A degeneracy in B, K and L (marginally 
also 5) becomes manifest even when looking at the radio spectrum as a function of 
the latitude. It is only when we turn to physical (i.e. B/C compatible) values of the 
propagation parameters that the degeneracy is partially broken in the spectrum w.r.t. 
the latitude, see Fig. [3j 

6. Conclusions and Outlook 

The connection between synchrotron radiation and radio data provides a very interesting 
means of inferring properties of the interstellar galactic electron population at low 
energies [33]. We have shown that the radio sky from MHz to GHz, when averaged over 
large scales, can be understood in terms of synchrotron emission of diffused galactic 
e predicted in a model consistent with many other cosmic observables, including CR 
nuclei (from protons to iron), radioactive isotopes and antiprotons. This is a rather new 
result and strong evidence that we are not too far from having a global picture of the 
phenomena occurring in the galaxy. 

We also found first indications that such a description breaks some of the 
degeneracies encountered when constraining the properties of the diffusive halo with 
nuclear CR data. A more detailed analysis of the allowed space of propagation 
parameters will be treated in a forthcoming publication [32]. However, let us stress 
that we have already presented preliminary evidence for lkpc < L < 15kpc. We note 
that in particular the p flux from dark matter annihilations in the galactic halo is mostly 
sensitive to the volume probed by the magnetic diffusion zone; a lower bound on L will 
thus have important implications for indirect dark matter searches. 
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